Extensive analytical and numerical investigation of the kinetic and stochastic Cantor 
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We investigate, both analytically and numerically, the kinetic and stochastic counterpart of the 
triadic Cantor set. The generator that divides an interval either into three equal pieces or into three 
pieces randomly and remove the middle third is applied to only one interval, picked with probability 
proportional to its size, at each generation step in the kinetic and stochastic Cantor set respectively. 
We show that the fractal dimension of the kinetic Cantor set coincides with that of its classical 
counterpart despite the apparent differences in the spatial distribution of the intervals. For the 
stochastic Cantor set, however, we find that the resulting set has fractal dimension df — 0.56155 
which is less than its classical value df = Nonetheless, in all three cases we show that the sum 
of the d/th power, df being the fractal dimension of the respective set, of all the intervals at all 
time is equal to one or the size of the initiator [0, 1] regardless of whether it is recursive, kinetic 
or stochastic Cantor set. Besides, we propose exact algorithms for both the variants which can 
capture the complete dynamics described by the rate equation used to solve the respective model 
analytically. The perfect agreement between our analytical and numerical simulation is a clear 
testament to that. 

PACS numbers: 61.43.Hv, 64.60.Ht, 68.03.Fg, 82.70Dd 



I. INTRODUCTION 



The history of describing natural objects by geome- 
try is as old as the history of science itself. Perhaps, 
the oldest of our pedagogical understanding about the 
properties of physical objects is geometry. Thanks to 
the Greek philosopher Euclid who is in fact the princi- 
pal architect for laying the early foundation of geometry 
which is now known as the Euclidean geometry. For cen- 
turies, it has been the only means of describing geome- 
try of physical objects. However, the then scientists also 
realized that the nature is not restricted to Euclidean 
space only; instead, most of the natural objects we see 
around us are so complex in shape that conventional Eu- 
clidean space is not sufficient to describe them. It was 
not until the work of Bcnoit B. Mandelbrot that dramatic 
progress was made. In 1975 Mandelbrot introduced the 
idea of fractal that has revolutionized the whole concept 
of geometry [l| . Prior to the inception of fractal, geom- 
etry remained one of the main branches of mathemat- 
ics. However, soon after its inception, it has attracted 
mathematicians, physicists, and engineers all alike and 
hence generated a widespread interest. All credit goes 
to Mandelbrot for the way he presented the idea of frac- 
tal through his monumental book The Fractal Geometry 
of nature Q. Indeed, he presented his book in an un- 
usually inspiring way and since then it remained as the 
most favourite standard reference book for both begin- 
ners and researchers. Due to its wide interest, it has 
brought many seemingly unrelated subjects under one 
umbrella and provided a tools to appreciate that there 
exists some kind of order even in the seemingly complex 
and apparently disordered many natural geometric struc- 



tures. 

The importance of fractal and multifractal in nonlinear 
dynamics especially in chaos can hardly be exaggerated. 
Their relation exist to such an extent that the impact of 
the ideas of chaos and fractals in physics and other scien- 
tific disciplines in recent years have been enormous. One 
common thing about both chaos and fractals is that much 
of their progress was possible due to high configuration 
computers and numerical simulations. Interestingly, the 
real progress of fractals and chaos both began during the 
1970s. Most of the books or articles written on chaos are 
found to invoke the idea of fractal since the plot of the 
fractal basins associated with chaos provides a qualita- 
tive idea of the extent of complications in the prediction 
of its future evolution Q . For this and because of their 
common history often general people regards the two con- 
cepts synonymous. Note that the keywords in chaos are 
non-linearity, unpredictability, sensitivity to initial con- 
ditions. On the other hand, the keywords in fractals are 
self-similarity and scale-invariance that is, the objects 
look the same on different scales of observation. 

The exact definition of fractal even after all these years 
is still elusive. Mandelbrot himself was somehow reluc- 
tant to confine it within the boundaries of a mere def- 
inition. He, nevertheless proposed that fractal can be 
defined as a geometric object which is similar to itself on 
all scales. That is, if one zooms in on a fractal object it 
will look exactly similar to the original shape. Mandel- 
brot offered also a mathematical definition "a fractal is 
by definition a set for which the Hausdorff-Besicovitch 
dimension strictly exceeds the topological dimension" . 
Sometimes fractal is also defined in terms of mass-length 
relation [4]. That is, if the mass M of an object is related 
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to its different possible size L by the relation 

M(L) ~ L df , (1) 

where df is less than the dimension of the space in which 
the object is embedded then the object is called fractal 
and it is quantified by its dimension df. 

The Hausdorff-Besicovitch dimension actually plays 
the pivotal role for a rigorous definition of fractal and 
for providing a procedure to find the fractal dimension. 
Consider that the measure Md is the size of the set of 
points in space of the object. We can quantify the size of 
the measure Md by using d-dimensional hypercube of lin- 
ear size 5 as a yardstick to find the number N(5) needed 
to cover Md 0). That is, we can cover the object to form 
the measure 

M d = Y,S d = N(6)5 d , (2) 

where S d is the test function. The number N(8) will 
obviously gets smaller as the size of the yardtick 6 gets 
larger. It can be easily shown that N(5) satisfies the 
following generalized relation 

N{S/n) = n dE N{5), (3) 

if Md describes an Euclidean object. The above relation 
for N satisfies the property of a homogeneous function 
N(Xx) = X p N(x) if one chooses A = n _1 , x = 6 and 
p = —ds- One can explicitly prove that only power-law 
solution can satisfy Eq. ([3]) 6]. Indeed, it is easy to 
check that N(S) ~ 5~ dE can solve Eq. J3J). In reality 
there exist another class of objects which also exhibits a 
similar power-law relation between the number N and the 
size of the yardstick S but with a non-integer exponent df 
instead of dE and hence one can generalize the relation 
as 

N{8) ~ S- D , (4) 

where D is called the Hausdorff-Besicovitch (H-B) di- 
mension Q. Using this relation in Eq. © we find that 
the value of Md assumes a finite constant value only if 
d = D otherwise the measure Md is either zero if d > D 
or infinity if d < D in the limit <5 — > 0. The H-B dimen- 
sion therefore is the critical dimension d = D for which 
the measure Md neither vanishes nor diverges as i5 — > 0. 
So, an object is fractal if the H-B dimension D = df a 
non-integer value. 

The best known text example of fractal is the triadic 
Cantor set. However, this classical Cantor set lacks in 
two ways from the natural fractals. Firstly, it does not 
appear through evolution in time although fractals in na- 
ture do so. Secondly, it is not governed by any sort of 
randomness throughout its construction. On the other 
hand, natural fractals always occur through some kind 
of evolution accompanied by some randomness. In this 
article, we therefore, investigate two interesting variants 
of the classical Cantor set, the kinetic and the stochastic 



Cantor set, in which time and randomness are incorpo- 
rated in a logical progression. The kinetic Cantor set 
differs from the classical Cantor set in the sense that at 
each step only one of the available interval is picked, with 
probability proportional to its size, for dividing it into 
three equal pieces and remove the middle third. On the 
other hand, in the stochastic Cantor set the definition 
of the generator is modified in the sense that it divides 
an interval into three pieces randomly and remove the 
middle third. However, the modified generator is applied 
sequentially to only one interval which is also picked ac- 
cording to its size at each generation step like in the 
kinetic Cantor set. The adventage of applying the gener- 
ator sequentially is that we can invoke time as one of the 
parameter in the problem. On the other hand modify- 
ing the generator in the stochastic Cantor set incorporate 
randomness into the problem as it evolves. In fact, it is 
well understood that in our world almost nothing is sta- 
tionary or strictly deterministic. Every natural objects 
are seemingly complex in character though they are gov- 
erned by simple rules since simple rules when repeated 
over and over again can appear mighty complex. 

The Cantor set provides us with a wealth of interest- 
ing properties which are taught in advanced undergradu- 
ate and graduate studies in discrete mathematics course. 
Apart from its pedagogical importance the Cantor set 
problem has also been of theoretical and practical inter- 
est. For instance, Sears et al has shown that a nonlinear 
system that supports solitons can be driven to generate 
Cantor set Q. In another case, it has been found that 
the electromagnetic wave is strongly enhanced and local- 
ized in the cavity of the Cantor set near the resonant 
frequency Q. Krapivsky and Redner shown that the 
probability distribution in a random walk with a shrink- 
ing steps is a Cantor set @ . Recently, Esaki et al stud- 
ied propagation of waves through Cantor set media and 
found some very interesting results such as complete re- 
flection or complete transmission including scaling prop- 
erty of the transmission co-efficients [To| ■ 

The remainder of this article is organized as follows. 
In section II, titled 'recursive Cantor set' we discussed 
the well-known triadic Cantor set. In section III, we pro- 
posed the kinetic counterpart of the triadic Cantor set 
and solved it exactly to obtain the fractal dimension and 
various other properties. We also proposed exact algo- 
rithm for the kinetic Cantor set to solve the model by 
numerical simulation. In section IV, we investigated the 
stochastic counterpart of the triadic Cantor set and once 
again we proposed its exact algorithm to solve it by nu- 
merical simulation. In section V, we summarized our 
results. 



II. RECURSIVE CANTOR SET (RCS) 

The novelty in the definition of the triadic Cantor set 
is its simplicity. The notion of fractal and its inherent 
character, self-similarity, is almost always introduced to 
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the beginner through this example. The triadic Cantor 
set can be defined as follows. It starts with an initia- 
tor of unit interval [0, 1]. The generator then divides it 
into three equal parts and deletes the middle third leav- 
ing behind two sub-intervals, each with size one-third of 
the original interval such as [0, ^] and [|, 1]. In the next 
step, the generator is again applied to each of the two 
sub-intervals that divides them into three equal parts of 
size g and remove the middle third from both. The pro- 
cess is then continued by applying the generator on the 
remaining intervals recursively ad infinitum and hence 
we call it recursive Cantor set (RCS). Like its definition, 
finding the fractal dimension of the RCS problem is also 
trivially simple. According to the construction of the 
RCS process there are N = 2" intervals in the nth gener- 
ation each of size S = 3~" and hence it is also the mean 
interval size. The most convenient yard-stick in the nth 
step, therefore, is the mean interval size S = 3~ n . The 
generation number n can be written as 

n = -InT (5) 

Using it in N — 2™ we find that the number N falls off 
following power-law against mean interval size S i.e., 

N(5)~5- d f, (6) 

with df = j^l . Since the exponent df of the above re- 
lation is non-integer and at the same time it is less than 
the dimension of the space d = 1 where the set is embed- 
ded, it is the fractal dimension of the resulting triadic 
Cantor set. Note that the set does not fill up the unit in- 
terval by uniform distribution of zero-dimensional points 
to describe a line rather, it fills up the unit interval by 
zero-dimensional points (also known as the Cantor dust) 
in such a special way that it possess exact self-similarity. 

III. KINETIC CANTOR SET (KCS) 

One may wonder what if we divide only one interval at 
each step instead of dividing every available intervals as 
done in the RCS problem? Clearly, the spatial distribu- 
tion of intervals along the line will be very different from 
the one created by the RCS problem. But, then the ques- 
tion is: Will the number N needed to cover the set by an 
yardstick, say of size S, still exhibit power-law against 81 
If yes, will the exponent vis-a-vis the fractal dimension 
be the same as for the RCS? To find a definite answer to 
these questions, the same generator that divides an inter- 
val into three equal pieces and remove the middle third in 
the RCS problem is applied to only one of the available 
intervals instead of applying it to each of the available 
intervals at each step. But after step one and beyond the 
system will have intervals of different sizes and hence it 
raises further question: How do we choose one interval 
when the system has more than one interval of different 
sizes? We choose the case whereby an interval is picked 



with probability proportional to their respective sizes as 
it appears to be the most generic case. One advantage of 
modifying the RCS problem in this way is that we can use 
the rate equation approach to solve it analytically. Since 
time becomes one parameter of the problem we call it 
kinetic Cantor set (KCS). 

The KCS problem can be defined as follows. Like RCS 
it may also starts with an initiator of unit interval [0, 1]. 
In the first step the generator therefore divides the ini- 
tiator into three sub-intervals of equal size and remove 
the middle third. The two newly created intervals are 
labelled as I and 2 starting from the left end of the two 
surviving intervals. In the next step we generate a ran- 
dom number R from the open interval (0, 1) and find 
which of the two subintervals I and 2 contains this num- 
ber R in order to ensure that intervals are picked accord- 
ing to their size. If R is found within the interval [0, -|] 
then we pick interval 1, if it is found within [|, 1] we pick 
interval 2. Say, the interval 1 contain R and hence we 
pick interval 1. The generator then divides it into three 
equal pieces and remove the middle third. The left end 
of the two newly created interval is then labelled with its 
parent label 1 and the intervals on its right is labelled 
as 3. In any case, time is increased by one unit even if 
R is found within the interval that has been removed. 
Note that between two successive generation steps the 
time unit may increase by several units since each time 
an attempt in picking an interval is unsuccessful the time 
is increased by one unit. One therefore cannot write a 
straightforward relation between time t and generation 
step j although we will explore later that there do ex- 
ist a non-trivial relation between N = (j + 1) and time 
t. The jth step therefore starts with j number of inter- 
vals labelled as 1,2, j and at the end of the jth step 
the system will have j + 1 intervals whose sizes can be 
denoted as x\, X2, Xj + \. The basic algorithm of the 
jth step which starts with j number of intervals can be 
described as follows. 

(a) Generate a random number R from the open inter- 
val (0,1). 

(b) Check which of the 1,2,..., j intervals contains the 
random number R. Say, the interval that contain R 
is labelled as m and hence pick the interval m. Else, 
if none of the j intervals contain R then increase 
time by one unit and go to step (a). 

(c) Apply the generator to the sub-interval m to divide 
it into three equal pieces and remove the middle 
third. 

(d) Label two newly created intervals starting from the 
left end which is labelled with its parents label m 
and the interval on the right is labelled with a new 
number j + 1 that has not been already used. 

(e) Increase time by one unit. 

(f) Repeat the steps (a)-(e) ad infinitum. 
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FIG. 1: Schematic illustration of the construction of the ki- 
netic Cantor set. On the left the numbers below j indicates 
generation steps which is related to number of intervals via 
N — j + 1. On the right the numbers below t indicates in- 
crease in time which is related to N via a non-trivial relation 
that we will explore. 



The expression for the mean interval size after the jth 
step therefore is S = J2i +± x i/(j + 1) an d the correspond- 
ing time t can be obtained from the counter used for time 
in the algorithm. 

To solve the KCS problem analytically we use the rate 
equation approach. The state of the system at any time 
t can be characterized fully by the interval size distribu- 
tion function C(x, t) which is defined in such a way that 
C(x, t)dx is the number of intervals in the size range x 
and x + dx at time t. The evolution of the distribution 
function C(x, t) can then be described by the following 
master equation 



dCjx.t) 
dt 



i-x-y 

-C(x,t) I dy / F(y,z\x)dz 
o Jo 

+ 2 / dyC(y,t)£ F(x,z\y)dz, (7) 



where the kernel F{x 1 y\z) describes the rate and rules 
how a parent interval of size z is divided into three smaller 
intervals of size x, y and (z — x — y). The first term on the 
right hand side of the above rate equation describes the 
loss of interval of size x due to breakup of an interval of 
size x while the second term describes the gain of interval 
of size x due to breakup of an interval of size y > x 
into three smaller pieces so that at least one of the three 
smaller intervals is of size x. The factor 2 in the gain 
term guarantees that only two of the three intervals are 
kept and one is removed which is exactly the definition 
of the Cantor set. Note that if the factor 2 in Eq. ([7|) 
is replaced by the factor 3, then the resulting equation 
describes the kinetics of ternary fragmentation process. 
The ternary fragmentation equation was first proposed 
and solved exactly by Ziff [12j]. In order to mimic the 
generator that picks an intervals according to size of the 
available intervals and divide it into three equal pieces 
we choose 

F(x, y\x + y + z) = (x + y + z)S(x - y)S(y - z). (8) 



where the two delta functions ensures that the three in- 
tervals produced by the generator are equal in size. Sub- 
stituting this kernel into Eq. ([7]) we get 

^^- = ~C(x,t) + 6xC(3x,t). (9) 

One can solve the above equation exactly to find the so- 
lution for C{x,t) using the Laplace transformation. 

We are not really interested in the solution for the 
distribution function C(x,t) rather we are interested in 
finding solution of its nth moment defined as 



M„(t) = / x n C(x,t)dx, 



(10) 



as it can provide almost all the information that we in- 
tend to find. Moreover, we find it more convenient to an- 
alyze the moment of the distribution function than the 
function itself. Incorporating the definition of the nth 
moment in Eq. @ yield the following rate equation for 
the nth moment 



Mn(t) 

dt 



3 n+l 



M n+1 (t). 



(11) 



It is interesting to note that the distribution function 
C(x,t) is not a directly observable quantity but its vari- 
ous moments are, for instance, Mo(t) = N(t) is the num- 
ber of available intervals at time t, Mi(t) = L is the sum 
of the sizes of all the available intervals at time t etc. This 
clearly justifies the reason behind focusing on finding the 
solution for the various moments than the distribution 
function itself. For consistency check, let us consider a 
case whereby none of intervals are removed after gener- 
ator divides an interval into three equal pieces. That is, 
the sum of the sizes of all the intervals at all time would 
be equal to the size of the initiator. The corresponding 
rate equation for M n (t) which can be obtained from Eq. 
(|11[) upon replacing the factor 2 by 3 from which one can 
easily find that 



dM x 
dt 



= 0, 



(12) 



and hence Mi (t) = L is indeed independent of time which 
is infact equal to the size of the initiator. 

In order to obtain a complete solution for the nth mo- 
ment of Eq. (fTT|) we assume that there exists a value 
n = n* so that M n * remains independent of time or 
a conserved quantity in the long time limit. Indeed, a 
closer look into the rate equation for the nth moment 
immediately implies that we can obtain the value of n* 
by applying the steady-state condition 



lim 



dM n (t) 



0. 



t — »oo dt 

in Eq. (fTT|) and hence obtain the following equation 

1 2 0. 



3 n* + l 



(13) 



(14) 
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Solving it for n* we immediately find that n* = 
which implies that M\n2 is a conserved quantity. One 
of the property of fractal is that it must obey scaling or 
self-similarity. As we are expecting that the KCS prob- 
lem like its cousin RCS problem will also generate fractal 
in the long time limit. It is therefore reasonable to an- 
ticipate that the solution of Eq. ([9]) for general n will 
exhibit scaling. Existence of scaling means that the var- 
ious moments of C(x,t) should have power-law relation 
with time and hence we can write a tentative solution of 
Eq. (JTT1) as below 



M n {t) ~ A(n)t a{n) . 



(15) 



If we insist that it must obey the conservation law, 

A/in2 = const, then we must have a(ln2/ In 3) = 0. Sub- 
in 3 \ i / 

stituting this in Eq. (jlip we obtain the following recur- 
sion relation 



a[n + 1) = a(n) — 1. 



(16) 



Iterating it subject to the condition that a (In 2/ In 3) = 
gives 



a(n) = -(n-— ). 



(17) 



We therefore now have an explicit asymptotic solution 
for the nth moment 



M n (t)~t 



where Z = — 1. 



(18) 



According to Eq. (fT8|) we find that the number of 
intervals N(t) — Mo(t) grows as 




FIG. 2: The decrease in mean interval size 5 against time t, as 
defined in the algorithm, is drawn in the log — log scale using 



numerical data and found that it satisfies 5 
also predicted by analytical solution. 



r 



which is 
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log(S) 



FIG. 3: Plots of log(iV) vs log(<5) is drawn using numerical 
data. The line has slope equal to df — revealing N ~ 8~ dl 
as predicted by the theory. 



N(t) ~t*». 



(19) 



and the mass or the sum of all the intervals size decreases 
against time as 



M(t) ~ t 



-0.369 



(20) 



The solutions for M (t) and M\(t) can provide us with 
the information how the mean interval size 5(t) — 
Mi(t)/M (t) decay and find that 



(21) 



where the kinetic exponent Z = — 1. We now apply the 
idea of fractal analysis which is briefly described in the 
introduction. To this end we find it convenient to use 
typical or mean interval size 8(t) as the yard-stick to mea- 
sure the resulting set since it will always give an integer 
N. This is equivalent to expressing the number N in 
terms of 5 and in doing so we find that N(5) decreases 
following the same power-law as Eq. ([6]) including its 
exponent df — The H-B dimension of the result- 

ing KCS problem therefore is D = which is exactly 
the same as its recursive counterpart albeit the spatial 
distribution is very different. 



To verify our analytical results we performed Monte 
Carlo simulation based on the algorithm (a)-(f). We 
first collected data for the mean intervals size 6 = 
X i/U + 1) against time where a^s are the size of 
the intervals specified by their labels i = 1, 2, This 
data is used to draw log(5) vs \og(t) in Fig. (2) and we 




FIG. 4: The sum of the d/th power of all the intervals x x + 

x 2 f + ...+Xj t — 1 regardless of time if df = |^-| which is equal 
to fractal dimension of the kinetic Cantor set. Its analytical 
counterpart is the d/th moment of the distribution function. 
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clearly find a straight line with slope equal to 1 reveal- 
ing that the mean interval size decreases following the 
same inverse power-law as predicted by Eq. (f2Tj) . Nu- 
merical data are also used to plot log(TV) against log(<5) 
and again we find a straight line but with a slope equal 
to which is exactly as predicted by our theory (see 
Fig. (3)). It clearly proves that the underlying mech- 
anism described by the Eq. ^ has been captured by 
the algorithm (a)-(f) we proposed for the KCS problem. 
One may think of another variant of the kinetic Cantor 
set by removing one of the three intervals randomly each 
time generator divides an interval into three equal inter- 
vals. Surprisingly, numerical data reveals that the value 
of df still is the same regardless of exactly which of the 
three intervals is removed each time an interval is divided 
into three equal intervals. Indeed, the rate equation for 
the distribution function C(x,t) given by Eq. ([9]) do not 
distinguish the three smaller intervals from one another. 
Now, incorporating df — in Eq. (jT5J) we can conclude 
that the moment of order equal to fractal dimension df is 
a conserved quantity. To verify this we collected data for 
the sum of the <i/th power of the size of all the intervals 
available at the generation step j and found 

xi f +x%'+xi' + +xf=l, (22) 

regardless of the value of j vis-a-vis time t which is equiv- 
alent to Md f ■ This is quite an extra-ordinary revelation. 
This result make us wonder if such conservation law also 
exist in the case of traditional recursive Cantor set. Due 
to recursive nature of the RCS problem, after any gener- 
ation step n all the N — 2™ intervals will have the same 
size 3~". The sum of the dfth power of all the intervals 
after the nth generation step therefore is 

2™(V n ) T ^ = l, (23) 

regardless of the value of n. We wonder what if we fur- 
ther modify the generator to create the stochastic fractal 
instead of random fractal. Will the relation that the sum 
of the dfth power of all the intervals at any given time 
be the same with that of at any other time? 

IV. STOCHASTIC CANTOR SET (SCS) 

The RCS problem hardly has any relevance to the frac- 
tal that appear in nature as it lacks at least in two ways 
from those that occur in nature. For instance, it does 
not have any kinetics but most natural fractal appears 
through some kind of evolution and it does not have any 
randomness but nature love to enjoy some degree of ran- 
domness to say the least. Though the KCS problem ap- 
pear through evolution but it still lacks in randomness. 
We therefore ask: What if we use a generator that di- 
vides an interval randomly into three smaller intervals 
instead of dividing into three equal intervals? To find 
an answer to this question consider an algorithm for the 
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FIG. 5: Schematic illustration of the construction of stochas- 
tic Cantor set. The numbers on the top of each intervals 
indicate their labels. The numbers on the left under j indi- 
cates generation steps and the numbers under t on the right 
illustrate how time evolves. 

stochastic Cantor set as described below. We start the 
process with an initiator of unit interval [0, 1] as before 
but unlike the previous cases the generator here divides 
an interval randomly into three pieces. The algorithm 
for the jth generation step that starts with j number of 
intervals can be described as follows. 

(i) Generate a random number R from the open inter- 
val (0,1). 

(ii) Check which of the 1,2, j intervals contains R. 
Increase time by one unit after every checking, 
starting from label 1 then label 2 and so on till 
an interval, say [a, b] labelled as k, is found. Go to 
step (i) if none of the available intervals contain R. 

(iii) Apply the generator onto the interval k to divide 
it randomly into three pieces. For this we generate 
two random numbers from the open interval (a, 6), 
say c and d where say c < d, to divide the inter- 
val into [a, c], [c,d] and [d, b] and delete the open 
interval (c, d) . 

(v) Update the logbook by labeling the left end of the 
two newly created interval [a, c] with its parents 
label k and right end of the two [d, b] is labelled as 
(J + 1). 

(vi) Increase time by two units since two cuts are needed 
to divide an interval into three smaller intervals. 

(vii) Repeat the steps (i)-(yi) ad infinitum. 

The generalized master equation for the Cantor set 
given by Eq. ([7]) can describe the rules of the SCS prob- 
lem stated in the algorithm (i)-(vii) if we choose 

F(x, y\x + y + z) = 1. (24) 

The master equation for the stochastic Cantor set then 
is 

^ = 4 C M + 2 J x ° C C(y,t)(y-x)dy. (25) 



log(S) 



FIG. 6: Decrease in the mean interval size 5 as a function 
of time t. The straight line in the log(5) vs \og(t) plot with 
slope equal to | implies that 5 ~ t~? and hence it is in perfect 
agreement with analytical prediction. 



This is exactly the rate equation first proposed and solved 
analytically by Krapivsky in the context of random car 
parking problem and later by an alternative method in 
the context of stochastic fractal [H, [l5| . In this article, 
we however give an exact algorithm that can capture the 
complete dynamics described by the above rate equation 
and verify the analytical results by numerical simulation. 
Following the method of Krapivsky and Ben-Nairn we 
incorporate the definition of the nth moment in the rate 
equation to give 



dM n (t) 
dt 



= 2 



and then solve it for M n (t). Following the same proce- 
dure as for the KCS problem we obtain the asymptotic 
solution for the nth moment 



with z = , 

2' 



(27) 



FIG. 7: Variation of the number of yard-sticks N needed to 
cover the stochastic Cantor set is shown against the size of 
the yard-stick 5. Data collected from simulation is averaged 
over 5000 independent realizations. Slope of this straight line 
is found to be equal to 0.56155 which means that N(S) ~ 

j-0. 56155 



FIG. 8: The sum of the d/th power of all the intervals x^ f + 

x 2 f + ... + Xj f = 1 regardless of time provided df = 0.56155 
which is the fractal dimension of the stochastic Cantor set. 
Data has been collected from the numerical simulation after 
averaging over 5000 independent realizations. 



where the number 0.56155 is the real positive root of the 
following quadratic equation 



, r(n + i)r(2) 

' r(n + 3) 



(28) 



Note that once again we find that the exponent of the 
power-law relation for M n (t) is linear in n and hence 
the system must obey a simple scaling but only in the 
statistical sense. It is interesting to note that the moment 
of order n — 0.56155, instead of n = In 2/ In 3 in the KCS 
problem, is a consserved quantity. Using Eq. (fT0|) in the 
definition of the mean interval size we find that 



5(t) ~t~*. 



(29) 



Once again we use it as the yard-stick and find that the 
number N(S) needed to cover the resulting set scales as 



N(S) ~ <r 



(30) 



where df = 0.56155 is the fractal dimension of the 
stochastic Cantor set. It reveals that the fractal dimen- 
sion of the stochastic Cantor set is less than that of its 
recursive or kinetic Cantor set. 

To verify our analytical results we performed numeri- 
cal simulation based on the algorithm (i) — (vii). Note 
that the definition of time in this algorithm is very much 
different from that of the kinetic Cantor set. Like in the 
KCS problem here too we collect data for the mean inter- 
val size S against time. We then plot log(J) versus log(i) 
in Fig. (6) and find straight line with slope equal ^. 
It implies that the mean interval size decreases following 
exactly in the same fashion as predicted by Eq. (|29|) . We 
also collected data for N against S. These data are used 
in Fig. (7) to show how the number N decreases with 
the yard-stick size S. A straight line in the logarithmic 
scale with slope equal to 0.56155 clearly implies that N 
exhibits power-law relation N ~ 5~ df with exponent df 
as predicted by Eq. (|3"0|) . Furthermore, in Fig. (8) we 
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show that the sum of the d/th power of all the intervals 

3,0.56155 + ^0.56155 + + + + ^0.56155 = j ^ 

regardless of the value of the generation step vis-a-vis the 
time. From the analytical point of view this is equiva- 
lent to the moment -M0.56155 of the distribution function 
C(x,t) which is indeed found to be independent of time. 
All these results clearly reveals that analytical results are 
in perfect agreement with the numerical simulation. 

V. SUMMERY 

To summarize, in this article we studied two interesting 
variant of the strictly self-similar triadic Cantor set. We 
solved the two models, the kinetic and stochastic Cantor 
set, both analytically using rate equation approach and 
numerically based on the exact algorithms that we pro- 
posed for both the problems. We found that the number 
of intervals and time are related via a generalized relation 
N ~ t df is true for both kinetic and stochastic Cantor set 
if we set d f = for the KCS problem d f = 0.56155 for 
the SCS problem. On the other hand, we found that the 



mean interval size decreases with time following power- 
law 5 ~ t z with exponent z = 1 for the KCS problem 
and Z = \ for the SCS problem. We then took the re- 
spective mean interval size as the yard-stick to measure 
the resulting set and found that the number TV needed 
to cover the set fall-off following power-law with expo- 
nent equal to their respective fractal dimension. We also 
found a generalized conservation law in the sense that 
the sum of the d/th power of all the available intervals at 
any time or the generation step is equal to the size of the 
initiator. This is true for recursive, kinetic and stochastic 
Cantor set. On the other hand, if we know the solution 
for the distribution function C(x,t), which can only be 
defined for kinetic and stochastic Cantor set, then the 
conservation law means that the d/th moment of C(x, t) 
remain independent of time. It is noteworthy to mention 
that such non-trivial conservation laws was also found re- 
cently in the context of condensation-driven aggregation 
process [l6j . We can perhaps conclude that emergence of 
fractal in a given system is always accompanied by some 
conservation laws which is ultimately responsible for fix- 
ing the various scaling exponents including the fractal 
dimension. 
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